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ABSTRACT 

We present observations of the peculiar supernova SN 1998bw, which was probably associated with 
GRB 980425. The photometric and spectroscopic evolution is monitored up to 500 days past explosion. 
We also present modeling based on spherically symmetric, massive progenitor models and very energetic 
explosions. The models allow line identification and clearly show the importance of mixing. From the 
late light curves we estimate that ~ 0.3 — 0.9 M© of ejected 56 Ni is required to power the supernova. 

Subject headings: supernovae: individual (SN 1998bw) — nucleosynthesis — gamma rays: bursts 



1. INTRODUCTION 

Supernova (SN) 1998bw was born famous because of its 
positional and temporal coincidence with the gamma-ray 
burst GRB 980425 (Galama et al. 1998a). The unique na- 
ture of the associated radio source, in terms of the high lu- 
minosity immediately after the explosion and the inferred 
relativistic expansion rate (Kulkarni et al. 1998), argues 
that SN 1998bw and GRB 980425 are related. Indepen- 
dent of the connection to the GRB, SN 1998bw is an inter- 
esting supernova. It was optically luminous and displayed 
very high expansion velocities. SN 1998bw was classified 
as a Type Ic (Patat & Piemonte 1998), and late time spec- 
tra reported here are consistent with this classification. 

Modeling of the early light curve and spectra suggested 
an extremely energetic explosion of a massive carbon- 
oxygen star (Iwamoto et al. 1998; Woosley, Eastman, 
& Schmidt 1999). As the energy was more than ten 
times that of a canonical core-collapse supernova, the 
term 'hypernova' was suggested for this event (Iwamoto 
et al. 1998). The mass of 56 Ni needed to power the early 
light curve in these models, 0.5 — 0.7 M , is much larger 
than the < 0.1 M©, typical for 'normal' core collapse SNe 
(Patat et al. 1994; Schmidt et al. 1994). This large mass 
of 56 Ni for SN 1998bw was disputed by Hoflich, Wheeler, 
& Wang (1999), who found that the early light curve can 
also be reproduced by an asymmetric SN Ic ejecting only 
0.2 M© of 56 Ni, given the right viewing angle and degree 
of asymmetry. 

The models for SN 1998bw differ considerably in terms 
of progenitor mass, nickel mass and explosion energy. We 
therefore focus on the late time evolution, when the nucle- 
osynthesis and density distribution can be studied directly. 
The total emission in the nebular phase is also less sensi- 
tive to asymmetries, compared to the early light curve. 

We present photometry and spectroscopy of SN 1998bw 
up to 500 days past explosion. We also compare these with 
results from spectral synthesis based on various progenitor 
models. A more detailed analysis will be given elsewhere. 



2. OBSERVATIONS, REDUCTIONS AND RESULTS 
2.1. Photometry 

Photometry was obtained at eight epochs between 33 
and 504 days past explosion, using the ESO-3.6m tele- 
scope on La Silla and the VLT/UT1 on Paranal. The 
images were bias subtracted and flat fielded using IRAF. 

To obtain the supernova magnitudes we used 
DAOPHOT. The instrumental magnitudes were converted 
using local standards from Galama et al. (1998b) 5 . At the 
earliest epochs, the errors are estimated as the standard 
deviations in the magnitudes of the local standards. They 
also encapsulate the neglect of color transformations. 

At later phases, the main uncertainty lies in the back- 
ground subtraction, as the supernova is superposed on an 
H II region. We constructed template PSFs using one or 
several stars, and the best PSF was chosen from judgments 
of the subtracted images. This PSF was subtracted from 
the SN position with a range of magnitudes, and the su- 
pernova magnitude was taken as the one which gave the 
smoothest background after subtraction. Upper and lower 
limits were estimated by noting when the subtracted im- 
age showed a hole, or a clear point source, at the supernova 
position. The supernova magnitudes are given in Table 1 
and plotted in Figure 1. 



2.2. Spectroscopy 

Spectroscopic observations were also obtained to late 
phases (Fig. 2). At ESO-3.6m we used EFOSC2 with 
grisms 11 and 12, and at VLT/UT1 we used FORS1 with 
the 300-grisms. All spectra were bias subtracted, flat 
fielded and wavelength calibrated. Flux calibration was 
done with a spectrophotometric standard star, and the 
absolute fluxing was performed using the simultaneous V- 
band photometry. All spectra were taken close to the par- 
allactic angle to reduce differential refraction. 

The spectral evolution of SN 1998bw from day 33 to day 
504 is shown in Figure 2. The evolution of the spectra in 
the nebular phase is fairly slow. The [O I] AA6300,6364 
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are the strongest lines, and the 7300 and 8600 A features, 
likely to be due to Ca II, decrease in strength relative to the 
[O I] lines. Other prominent features are seen at <~ 4570 
and 5890 A, and are probably due to Mg I] (or [Fc III]) 
and Na I D (and possibly He I), respectively. 

3. DISCUSSION 
3.1. General considerations 

A lower limit of 0.22 M of 56 Ni for SN 1998bw was 
achieved by McKenzie & Schaefer (1999) by comparing the 
^-magnitude at 170 days to that of SN 1987A. While the 
light curve of SN 1987A closely followed the decay rate of 
56 Co, suggesting full trapping of the 7-rays, the light curve 
of SN 1998bw fell substantially faster. The assumption of 
full trapping at 170 days in SN 1998bw used by McKenzie 
& Schaefer is clearly too conservative. In fact, the light 
curve of SN 1998bw is similar to that of many other SNc 
Ib/c (Clocchiatti & Wheeler 1997; Sollerman, Leibundgut, 
& Spyromilio 1998), and can be reproduced by a simple 
model of a radioactively powered ejecta leaking 7-rays due 
to homologous expansion (Fig. 1). If the decay would fol- 
low full trapping from day 64 and onward, the flux in the 
V-band at day 170 would be 2.4 times brighter than actu- 
ally observed (Fig. 1). Assuming full trapping at day 64, 
at the beginning of the fast decay, then provides an upper 
limit of ~ 41% on the trapping at day 170, increasing the 
lower limit of the Ni-mass to ~ 0.5 M©. 

The above estimate depends, however, on the bolometric 
correction and its evolution for these supernovae. That the 
bolometric decay rate is similar to that of the y-band at 
this epoch is indicated by the simple estimates of the bolo- 
metric luminosity by McKenzie & Schaefer (1999). This is 
also confirmed by our detailed modeling (see below) and 
supported by the fit to the y-curve by the simple 7-leakage 
model. However, the bolometric correction could still be 
different for SNe 1987A and 1998bw. If SN 1998bw emit- 
ted a smaller fraction of its bolometric luminosity in the 
optical range than SN 1987A, its Ni-mass would be over- 
estimated. 

Adopting a distance of 35 Mpc (Ho = 72), a reddening 
of E(B — V) = 0.06 (Schlegel, Finkbeiner, & Davis 1998), 
and integrating our spectrum at day 141, gives a lumi- 
nosity (L op t = 3.5 x 10 41 ergs s _1 ) that requires 0.1 M Q 
of 56 Ni. This limit clearly has to be corrected upwards 
as it assumes full 7-ray trapping, and that our spectrum 
contains all the flux emitted by the supernova. With the 
same argument as above, no more than 48% of the en- 
ergy is deposited at 141 days. Furthermore, detailed mod- 
els (see below) show that the observed spectrum contains 
only 71% of the total flux. These corrections indicate that 
the lower limit must be increased to 0.3 M©. The domi- 
nant errors are the uncertainties in bolometric correction 
(~ 10%), distance, and reddening for SN 1998bw. Allow- 
ing for a 15% uncertainty in the distance scale, and ±0.1 
mag in A v , make the absolute fluxes (and the Ni-mass) 
uncertain by ~ 35%. 

3.2. Detailed modeling 

To obtain more quantitative estimates of the Ni-mass, 
detailed modeling is required. Our analysis uses an up- 
dated version of the spectral code in Kozma & Fransson 



(1998a, 1998b). At the epochs discussed here, the ejecta 
emission is powered by the radioactive decay of 56 Co. The 
deposition of the resulting 7-rays and positrons is mod- 
eled in detail according to Kozma & Fransson (1992). The 
ionization and heating is balanced by thermal and radia- 
tive processes, producing the observed optical and infrared 
emission. This time-dependent code gives good agree- 
ment with light curves and late time spectra for SN 1987A 
(Kozma & Fransson 1998a, 1998b; Kozma 2000). 

As input models we use different combinations of un- 
mixed and artificially mixed models of SN 1998bw, based 
on the C0138-model calculated by Iwamoto et al. (1998) 
and the C06-modcl from Woosley et al. (1999). Both 
models are spherically symmetric explosions of massive 
carbon-oxygen stars. The C0138-model has a mass Mco 
= 13.8 M©, explosion energy i? cxp = 3 x 10 52 ergs, and 
mass of 56 Ni, M( 56 Ni) = 0.7 M Q . Such a carbon-oxygen 
core corresponds to a zero age main sequence (ZAMS) 
star of M Z ams ~ 40 M Q . The C06C-model with M co = 
6.55 M©, £ cxp = 2.2 x 10 52 ergs, and M( 56 Ni) = 0.47 M© 
corresponds to a Mzams ~ 25 M© star. 

In Figure 3 we show our modeled spectra together with 
the observations at 141 days. From the observed line 
profiles we conclude that mixing is important. Calcula- 
tions based on non-mixed models give very broad and flat- 
topped line profiles which are not seen in the observations. 
In the mixed models, we have also decreased the expansion 
velocities to get a better agreement with the observed line 
widths. A mixed model with the original velocities results 
in nicely peaked, but too broad, line profiles. A decrease in 
velocity results in an increase of the densities given by the 
original explosion models. We have increased the densities 
by factors of 10 - 100, except in the Fe-rich regions where 
it is kept unchanged (see below). These changes also affect 
the kinetic energy in the modified models. In the layers 
below ~ 12 000 km s _1 , which is the only material which 
contributes significantly to the emission at late times, the 
kinetic energies are 2 x 10 51 ergs (C06) and 5 x 10 51 ergs 
(C0138), respectively. Note that our analysis is not sen- 
sitive to the high velocity gas, which might also contain a 
significant amount of energy. Early time analysis is there- 
fore more reliable for determining the total kinetic energy 
of the explosion. As seen in Figure 3, the synthetic line 
profiles are still somewhat broader than the observed. We 
can overcome this by decreasing the velocities even more, 
but the high density in the oxygen-rich regions then re- 
sults in strong O I lines at AA 5577, 7774, 9265, which are 
not seen in the observations. 

Most features in the observed spectrum have a counter- 
part in the modeled spectra, and in Figure 3 we include 
some of the line identifications. The 4000 — 5500 A regime 
is problematic, and preliminary results from improved 
models suggest that most of these features are a mixture of 
Fe II and Fe III lines. The models presented here have too 
low a degree of ionization, with mostly Fe I and Fe II emis- 
sion. A higher density in the Fe-emitting material would 
result in even lower ionization, which is why we do not 
increase the density in the Fe-rich regions. The ionization 
can instead be increased by decreasing the density in these 
regions. Scattering and fluorescence may also play a sig- 
nificant role in understanding these features. We therefore 
regard the Fe I emission as dubious. 
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3.3. The Light Curve 

In Figure 1 we show the light curves from our calcula- 
tions of the mixed C06- and C0138-models (solid lines), 
and from the original, non-mixed models (dotted lines). 
The C0138-model is brighter than C06. This is not 
only due to the higher Ni-mass, but also a result of the 
higher ejecta-mass (M c j), and lower expansion velocity 
(u max ), giving a larger optical depth to the 7-rays (r 7 ); 
t 7 oc M cj /^ ax . 

The effect of changing the expansion velocities can also 
be seen in Figure 1. While the two C06-models have 
the same Ni-mass and ejecta-mass, they differ in expan- 
sion velocity. The mixed models are somewhat brighter 
than the non-mixed models, due to the higher 7-ray op- 
tical depth for the models with lower velocities. In the 
original C06- model we find that at 200 (400, 600) days 
only ~ 7.2 (4.4, 3.9)% of the energy released in the 
56 Co-decays is deposited in the ejecta. Out of this en- 
ergy ~ 47 (78, 87)% is from the kinetic energy of the 
positrons. For the mixed C06-model with decreased ve- 
locities, ~ 11 (5.7, 4.5)% of the energy is deposited. 

From the light curves, and within the framework of 
our modeling, we conclude that the nickel mass needed 
to power the late emission is indeed very high. For the 
two models that best reproduce the late time observations, 
we calculate the deposition of the radioactive energy, and 
thus the required mass of nickel. The C06-model requires 
0.9 M Q of 56 Ni to maintain the luminosity at 340 - 410 
days. The C0138-modcls, which traps 10% of the energy 
at 400 days, requires 0.5 M . The errors on these esti- 
mates are dominated by the uncertainty in distance and 
reddening stated in §3.1. 

Hoflich ct al. (1999) found that the early light curve of 
SN 1998bw could be explained in terms of a more 'normal' 
SN Ic, which ejected 0.2 M Q of 56 Ni, given the right degree 
of asymmetry. However, at the nebular phase studied here 
the ejecta are optically thin, and for a given density struc- 
ture the luminosity is not as sensitive to asymmetries as is 
the early light curve. By modeling the late time line pro- 
files, providing a measure of the energy input with velocity, 
we get a handle on the distribution of the ejecta, and can 
therefore model the gamma-deposition self-consistently. 

The Ni-mass determined for the C0138-modcl is close 
to the lower limit from §3.1. This is because this massive 
model closely represents the maximum trapping allowed by 
the observed light curve. Also, as the mixed C06-model is 
already dominated by positrons at 400 days, the estimated 
Ni-mass for this model must be close to the upper limit. 

While SN 1998bw ejected a large mass of 56 Ni, a few 
other supernovae have shown very low amounts of ejected 
56 Ni (SN 1994W 0.015 M , Sollcrman, Cumming, & 
Lundqvist 1998] and SN 1997D [~ 0.002 M Q , Turatto et 
al. 1998]). The iron yield from core-collapse supernovae 
thus seems to vary substantially. 

At very late phases (> 400 days), the low energy input 
means low temperatures for the ejecta, and most of the 
emission in our models comes out in the infrared. The op- 
tical light curves no longer directly follow the bolometric 
light curve. After 400 days the model light curves drop 
faster than the observations. A reduced density in the 
Fe-rich material (§3.2) can increase the temperature and 
boost in particular the F-curve at late phases. 



Another possible source for the powering at late phases 
is interaction with a circumstellar medium (CSM). The 
absence of H and He in the spectra means that substantial 
amounts of such gas should be present. Models of the radio 
emission indicate that some interaction is already taking 
place (Li & Chevalier 1999), although the CSM density in 
these models is probably too low to affect the optical light 
curve, as the wind velocity is high (~ 2000 km s^ 1 ). We 
see no significant spectroscopic signature that indicates 
CSM interaction, in our late spectra. 

4. SUMMARY 

We present photometry and spectroscopy of SN 1998bw 
up to 500 days past explosion. Starting with two differ- 
ent spherically symmetric explosion models, we calculate 
late time spectra and light curves using a detailed spectral 
synthesis code to model our observations. We conclude 
that mixing is important, and we also have to alter the 
density distribution macroscopically to obtain reasonable 
agreement with the line profiles in the observed spectra. 

Representing the supernova density profile with the dis- 
tribution which best reproduces the late time observations, 
we calculate the deposition of the radioactive energy, and 
determine how much 56 Ni is required to produce the late 
time luminosity. We find that a large amount, roughly half 
a solar mass, is needed to power the light curve one year 
after the explosion, although the exact number depends 
on the chosen explosion model. 

Evidence for asymmetries in SN 1998bw comes from 
early polarization measurements (Kay et al. 1998). The 
modeling presented here is restricted to the available sym- 
metric models, although the importance of mixing may 
point to an asymmetric explosion. Whether asymmetric 
models can reproduce the observations is of interest for 
future studies. 
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Table 1 
supernova magnitudes 



Epoch a Telescope V R I 



33 


3.6m 


14.69±0.05 


14.31±0.05 


14.26±0.05 


141 


3.6m 


17.34±0.05 


17.06±0.05 


16.78±0.05 


215 


3.6m 


18.63±0.05 


18.10±0.05 


17.92±0.05 


340 


3.6m 


20.7t^ 


19.9±0.15 


19.65±0.2 


357 


VLT 


21 +0 - 25 

Zi - U -0.15 


20 2+ 15 


20.0±0.15 


414 


VLT 


21.7±0.3 


21.0±0.20 


20.9+0^ 


471 


3.6m 


21 q +0 - 5 


21.6±0.30 


21 5 +o:3 


504 


VLT 


99 o+0. 4 
ZZl J -0.2 


22.1±0.30 


21.8±°; 5 2 



a Epoch in days past April 25, 1998, which we take as the date 
of explosion. 
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Days past explosion Days past explosion 



Fig. 1. — The upper left panel shows a simple model (solid line) compared to the F-band data described below. The flux 
is given by e -t / m - 3 x (1 —0.965 e~ r ), where 111.3 days is the decay time of 56 Co, and the 7-ray optical depth r = (t\/t) 2 
(here t\ = 120 days). The dashed curve is the 7-ray contribution and the dotted line the positron contribution. The 
dashed-dotted line shows the expected flux for full trapping from day 64 and onward. The other panels show the V, R, 
and / light curves, respectively, including, in addition to our data (large empty circles), observations from Galama ct al. 
(1998a) (1- 57 days past explosion), McKenzie & Schaefer (1999) (days 63 - 187), and Patat et al. (2000) (days 324 - 
417). The models included are C0138 and C06 with mixed ejecta and decreased velocities (solid lines), and the original 
C0138- and C06-models (dotted lines). The Am = ±0.39 error bar in the F-panel shows the uncertainty in distance and 
reddening. The thick solid line in the i?-panel shows the decay rate of 56 Co for full 7-ray trapping. 
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Fig. 2. — The spectral evolution of SN 1998bw from day 33 to day 504. For clarity the spectra have been shifted by 
a constant factor, A log F\ : +33 d (0.0), 141 d (-0.15),' 215 d (-0.5), 414 d (-1.0), 471 d (-2.0), and 504 d (-3.2). The 
wavelength scale has been corrected for the velocity of the parent galaxy, 2550 km s _1 . 
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Fig. 3. — Synthetic spectra based on the C0138-model (upper) and the C06-model (middle) together with the observed 
spectrum (lower) at 141 days. The C0138- and C06-models were artificially mixed with decreased velocities as described 
in the text. For some of the line identifications we indicate in which of the two models the line is important. The Fe I 
and Fe II lines are strongest in the C0138-model. The synthetic spectra have been shifted vertically by A log F\ = -1.5 
(C06) and A log F\ = -2.0 (C0138). 



